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1 . FOREWORD 


The Solar Energy System Economic Evaluation - Final Report has been 
developed by the George C. Marshall Space Flight Center as a part of 
the Solar Heating and Cooling Development Program funded by the 
Department of Energy. The analysis contained in this document de- 
scribes the economic performance of an Operational Test Site (OTS). 

The objective of the analysis Is to report the long-term economic per- 
formance of the system at Its installation site and to extrapolate to 
four additional locations which have been selected to demonstrate the 
viability of the design over a broad range of environmental and economic 
conditions. 

The contents of this document are divided into the following topics: 

• System Description 

• Study Approach 

• Economic Analysis and System Optimization 

• Results of Analysis: Technical and Economic 

• Economic Uncertainty Analysis 

• Summary and Conclusions 

The data used for the economic analysis have been generated through eval- 
uation of the Operational Test Site described in this document. The data 
that have been collected, processed, and maintained under the OTS Develop- 
ment Program provide the resource from which inputs to the simulation 
programs used to perform technical and economic analysis are extracted. 

The Final Report document, in conjunction with the Seasonal Report [3]* for 
each Operational Test Site in the Development Program, culminates the technical 


♦Numbers in brackets designate references found in Section 8. 
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activities which began with site selection and instrumentation system 
design in April, 1976. The Seasonal Report emphasizes the technical 
analysis of solar systems performance. It compares actual performance 
with predicted performance derived through simulation methods where 
actual weather and loads defined the inputs. The simulation used for 
final report analysis is based on the technical results of the seasonal 
report simulation, with the exception that long-term weather, and de- 
rived loads are used as inputs instead of measured weather and loads. 
This causes the expected value of solar system performance in the 
Seasonal and Final Reports to differ. In addition localized and stan- 
dard economic parameters are used for economic analysis in the final 
report evaluation. The details of the simulation program are described 
In References [4] and [5], Other documents specifically related to the 
solar energy system analysed in this report are [1] and [2], 
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2. SYSTEM DESCRIPTION 


The Solaron Akron Solar Energy System was designed to provide both space 
heating and domestic hot water (DHW) preheating for a dual level single- 
family residence containing approximately 1840 square feet in Akron, Ohio. 

Solar energy collection is accomplished with flat-plate collectors using 

air as the transport fluid. The collector array has a gross area of 546 
square feet and faces south at an angle of 45 degrees from the horizontal. 

Solar energy is stored in a 270 cubic foot rock thermal storage bin located 
on the lower level of the house. Solar energy is transferred to the DHW 
subsystem by means of an in-duct heat exchanger (HX1) whenever the system 
is storing collected solar energy. Water from the 80 gallon preheat tank 
and make-up water are transferred from the preheat system to the 52 gallon 
DHW tank when there is a demand for hot water. The auxiliary space heating/ 
cooling subsystem consists of an air to liquid heat pump coupled with a 1000 
gallon water storage tank. The heat pump can provide energy either directly 
to/from the house or the 1000 gallon tank in the appropriate season. The 
system is designed so that the heat pump can charge the 1000 gallon tank during 
off-peak hours when electrical rates are lower. Energy stored in the tank can 

then be used for space heating purposes as required. Chilled water stored in 
the tank can be used for space cooling as required. Auxiliary energy for both 
the space heating and DHW subsystems is provided by electricity. The heat pump 
has a nominal capacity of 30,000 Btu/Hr with supplemental heat strips rated at 
12 kW, and the auxiliary hot water heater is rated at 4.5 kW. The system is 
shown schematically in Figure 2-1, and sensor designations in Figure 2-1 are in 
accordance with NBSIR-76-1137 [6]. The measurement symbol prefixes: W, T, 

EP, and I represent respectively: flow rate, temperature, electric power, 

and insolation. The system has the following modes of operation: 

A. First Stage 

1. Collector to Storage and DHW. In this mode the collector blower 
transfers solar energy from the collector array to the rock 
thermal storage bin through the DHW heat exchanger. Part of the 
solar energy is utilized in the DHW preheat loop and the remain- 
ing solar energy is delivered to storage. This mode is entered 
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Figure 2-1 Solaron-Akron Solar Energy System Schematic 







whenever the differential temperature between the collectors 
and the return air duct is 40 + 7°F and heating demands are 
such that direct space heating from the collector array is not 
required. This mode terminates whenever the differential tem- 
perature falls to 25 + 5°F, or less, or direct space heating 
from the collector array is required. 

2. Collector to Space Heating Load. In this mode dampers MD1 and 
MD2 are open and solar energy goes directly to the residential 
area utilizing both the collector and circulating blowers. The 
DHW heat exchanger is bypassed in this mode and all collected 
energy is delivered to the space heating load. The same differ- 
ential temperature conditions described above also control opera- 
tion in this mode. 

3. Storage to Load. When Incident solar energy on the collector 
array is insufficient, space heating is provided from the 
storage bin by way of the circulating blower. Dampers MD1 and 
MD3 are closed in this mode and MD2 is open. A minimum storage 
temperature of 90°F is required for operation in this mode. 

B. Second Stage 

4. Heat Pump Auxiliary Direct. When insufficient solar energy is 
present on the collector array and the storage temperature is 
also insufficient to maintain a level of comfort, dampers MD1 
and MD2 close and MD3 opens to provide heated air from the heat 
pump by way of the auxiliary heating/cooling heat exchanger. 

At outdoor temperatures of approximately 40°F or above, the 
heat pump will carry the entire space heating load. For tempera- 
tures between 2°F and approximately 40°F, the heat pump is supple- 
mented by the electrical strip heaters. 

It is also possible to heat in this mode while, at the same time, 
collected solar energy is being delivered to storage. This 
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condition exists whenever the room thermostat is calling for 
second stage heating and sufficient insolation is available 
to allow the collector array to operate. 

5. Auxiliary Heat from Heat Pump Storage. This mode allows space 
heating from the off-peak water storage tank. During off-peak 
hours, when the heat pump is not needed to heat the residence, 

it stores hot water for use during this mode. Dampers MD1 and 
MD2 are closed and MD3 is open in this mode. 

C. Third Stage 

6. Electrical resistance (strip) heat is used whenever the heat 
pump is unable to maintain the desired comfort level in the 
house. Above 2°F the strips supplement the heat pump, as 
described in Mode 4 above, and below 2°F the strips carry the 
entire load. 

D. Space Cooling 

7. This mode in unrelated to solar energy system operation, but 
provides a unique capability to the auxiliary system which 

is intended to reduce operating expenses. When the heat pump 
is operating in the air conditioning cycle, heat energy can 
be removed from the house directly or from the 1000 gallon 
storage tank during hours when off-peak electrical rates are 
in effect. The chilled water in the storage tank is then 
pumped through the auxiliary heat exchanger (HX2) to supply 
cooling in the place of direct operation of the heat pump 
operating in the air conditioning cycle during peak hours. 


6 



3. 


STUDY APPROACH 


3.1 Introduction 


The Final Report is an economic evaluation of the solar energy system 
Chased on life cycle costs versus energy savings) for five cities which 
are considered to be representative of a broad range of environmental and 
economic conditions in the United States. Life cycle costs provide a mea- 
sure of the total costs of owning and operating a system over the life of 
the system rather than focusing solely on the initial cost of the system. 

The life cycle costs used in this evaluation consider hardware, instal- 
lation, maintenance, and operating costs for the solar-unique components 
of the total system. Energy savings result from replacement of conven- 
tional forms of energy by solar energy after the costs of producing the 
solar energy are deducted. The total system operates in a scenario that 
comprises long-term average environmental conditions, loads, fuel costs 
and other economic factors that are applicable in each of five cities. 

The five cities include four standard analysis sites which were selected 
according to the criteria listed below and the site where the system was, 
in fact, installed and operated. The selection criteria were based on: 

0 Availability of long-term weather data 
t Heating degree days (load related factor) 

0 Cold water supply temperature (load related factor) 

0 Solar insolation 

0 Utility rates 

0 Market potential 

0 Type of solar system 

To achieve the range of environmental and economic parameters desired, 
the four locations listed below plus the actual installation location, 
were used. A solar energy system buyer may evaluate his own local environ- 
mental and economic conditions relative to those considered in this Final 
Report by comparing the insolation available, the heat load, and the utility 
rates against the results reported in Section 5. 
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Albuquerque, NM 


2 

1828 Btu/Ft -Day average insolation* 

Medium heating load (4292 Heating Degree Days (HDD)) 

High utility rates (> 0.06 $/kWh )** 

Fort Worth, TX 

2 

1475 Btu/Ft -Day average insolation* 

Light heating load (2382 HDD) 

Medium utility rates (0.04-0.06 $/kWh )** 

Madison, WI 

2 

1191 Btu/Ft -Day average insolation* 

High heating load (7730 HDD) 

Medium utility rates (0.04-0.06 $/kWh )** 

Washington, DC 

2 

1208 Btu/Ft -Day average insolation* 

Medium heating load (5010 HDD) 

High utility rates (> 0.06 $/kWh)** 

Akron, Ohio 

2 

1086 Btu/Ft -Day average insolation* 

High heating load (6224 HDD) 

Medium utility rates (0.04-0.06 $/ kWh ) ** 

The parameters that define the system design were derived from the actual 
operating conditions of the system at the installation site. Solar energy 
system design may be economically optimized for the site at which the 

*Insolation values are average daily long-term values on a horizontal 
surface. 

**Utility rates are effective year-round averages based on 1000 kWh for 
Jan. 1980. See Appendix D. 
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system Is Installed. The fundamental objective in optimizing the design 
of a solar energy system on an economic basis is to minimize cost by 
allocating the required amount of energy between the solar and conventional 
portions of the system. To attain this objective, each unit of energy 
should be produced by the portion of the total system which generates 
the lowest Incremental cost In producing that additional unit of energy. 

This is accomplished In the final report analysis by determining the 
optimal solar energy system size (collector area or equivalently, solar 
fraction). 

In the Operational Test Site (OTS) Development Program there are many solar 
energy systems designed by many different contractors. Some of the designs 
were Installed in new buildings and some were retrofitted to existing build- 
ings. Consequently, there are a variety of factors which contributed to the 
design of a system at a given site. In some cases the objective of optimizing 
the design according to the previously stated criterion could not be met. A 
method of evaluation which establishes a common basis for evaluation of all 
these systems was required. The method selected is to optimize the collector 
size through the f-Chart [4], £5] design procedure. F-Chart is a design 
program developed by the University of Wisconsin for solar heating and/or do- 
mestic hot water systems. The program uses a set of design charts (developed 
by detailed simulations) which estimate the thermal performance of a solar 
system based on collector characteristics, storage, energy demands, and 
regional long-term weather data. Using the results of thermal analysis, an 
Iterative procedure Is implemented to select a collector area which minimizes 
the life cycle costs. Once the optimal collector size has been determined, 
the resulting thermal and economic performance can be obtained. 

The resolution of two inter- related problems was required In order to adapt 
f-Chart to the evaluation developed in the Final Report. The first was how 
to use the data and experience gained from the actual operation of the solar 
energy system; the second was what procedure to follow in view of the fact 
that all solar energy systems to be analysed do not have optimal collector 
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area sizing. To resolve the first problem, the characteristics of 
design and operation of the existing solar energy system were used to 
develop the input parameters for f-Chart. This procedure, detailed in 
Appendix A, involved the normalization of collector flow rates and stor- 
age capacity to collector area. Collector characteristics developed 
from field data through a collector analysis program were substituted for 
the theoretical single panel parameters furnished by the collector manu- 
facturer. To resolve the problem of different collector areas, an optimal 
collector area was derived for each site. The final adaption of f-Chart 
includes the inputs derived from operational data and optimal collector area. 

In addition to the f-Chart problems described above, certain internal 
modifications were required to enable the economic analysis of space 
heating and domestic hot water systems where the auxiliary energy source 
for space heating was a heat pump. This involved modification of the 
loads from which the economic parameters were computed. To modify the 
loads two coefficients of performance, i.e., SHCOP for the space heating 
system and HWCOP for the hot water system, which are described in Appendix 
A, were introduced. These COP's are used to adjust the cost of auxiliary 
energy considering the efficiency of the hot water system (assumed to be 
100 percent, therefore HWCOP =1), and the space heating system with its 
heat pump auxiliary. (See Table 5.1-3 for SHCOP at each analysis site). 

As the system application at each of the five analysis sites is studied, 
the loads are iteratively redefined, the site peculiar parameters are 
changed as described in Appendix A, and a new optimal collector area is 
computed. The economic factors are the result of the f-Chart analysis 
with these inputs. 
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3.2 Groundrules and Assumptions 


The cost differential between solar and the conventional system is 
significant to the economic evaluation in the Final Report. Cost 
items which were equal for both alternatives do not contribute to the 
differential cost. The cost of the conventional system was assumed 
to be identical with or without the solar alternative. Although a con- 
ventional system is usually selected according to the availability and 
cost of energy In a particular geographic region, this alternative is 
not permitted In the final report analysis because an existing system 
Is being evaluated. Savings which might be realized by comparing solar 
against an auxiliary other than the design option were not evaluated. 

The system configuration, including the conventional auxiliary, is the 
same for all five analysis sites. 

The cost of the solar-unique hardware is based on mass production esti- 
mates. The total incremental costs for acquisition of a solar alterna- 
tive are the sum of a cost proportional to collector area and a cost 
independent of collector area. For economic evaluation, life cycle 
costs (i.e., costs of acquiring, operating and maintaining the solar 
systems) were forecast on an annual basis over the design lifetime of 
the system, then discounted to an equivalent single constant dollar 
(1980) value as described in Section 4. 

Fuel costs are calculated at current (1980) local values for each of the 
five analysis sites. Other economic parameters are standardized by 
referencing current national economic conditions. Maintenance, insurance, 
depreciation, system life, salvage values (for commercial systems) are 
determined from best experience. Tax credits allowed by the Federal 
Government for the solar energy systems are credited against the acquisi- 
tion cost. A combined state and federal Income tax rate of 30 percent 
is assumed for estimating tax savings resulting from the interest paid 
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In financing a solar system. Property taxes arising from the increased 
value of property with an installed solar system are neglected due to 
the current trend In many states to forego these taxes to prevent them 
from being a disincentive to solar energy usage. 

The primary measure of cost effectiveness of the solar system for the 
evaluation In the Final Report is: 

t Life Cycle Cumulative Savings (LCCS) - The present value of the 
cumulative energy savings (in dollars) that result from operation 
of the solar system instead of the conventional system. 

Two secondary measures that depend on life cycle cumulative savings 
are: 

t Year of Positive Savings - Year in which solar system first 
becomes profitable; i.e. , the annual conventional fuel bill 
without solar exceeds the sum of the annual fuel bill with 
solar and the annual cost for the solar system. 

e Year of Payback - Year in which the compounded net savings 
equals the initial cost for the solar system. Net savings 
are computed with respect to the fuel cost of the conven- 
tional system. 
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4. ECONOMIC ANALYSIS 


4.1 Factors in Life Cycle Costs and Savings 

The economic calculations of this study are performed in the f-Chart 
program and are based on comparisons of life cycle costs of conventional 
energy systems with those of solar energy systems. The life cycle sav- 
ings of a solar energy system over a conventional energy system can be 
expressed as the difference between the total fuel savings that result 
from operation of the solar energy system and the increased costs that 
result from the investment in, the operation of, and maintenance of the 
solar energy system. The savings can be expressed by the relationship [7] 

LCCS = P 1 (C fe L e + C fe L f /C 0P f )F - P 2 (C A A + C £ ) (1) 

where LCCS = Life cycle cost savings of the solar 

energy system ($) in terms of present worth 

P-j = Factor relating life cycle fuel cost savings 
to first year cost savings 

C FE = Electrical energy cost per unit ( $/Mi 1 1 ion Btu) 

COPp = Heating system coefficient of performance 
L E = Hot water load (Million Btu) 

Lp = Space heating load (Million Btu) 

F = Solar fraction 

?2 ~ Factor relating life cycle investment 
operation and maintenance expenditures 
to the initial investment 

C. = Solar energy system costs dependent 
M 2 

on the collector area ($/ Ft ) 

A = Collector area (Ft^) 

C E = Solar energy system costs that are independent 
of collector area. ($) 
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It was assumed that the costs of components which were common to both 
conventional and solar heating systems (e.g. the furnace, ductwork, 
blowers, thermostat), and the maintenance costs of this equipment, are 
Identical. Consequently, all references to solar energy system costs 
refer to the cost increment above the common costs. 

The multiplying factors, PI and P2, facilitate the use of life cycle 
cost methods In a compact form. Any cost which was proportional to either 
the first year fuel cost or the initial investment can be included. These 
factors allow for variation of annual expenses with inflation and reflect 
the time value* of money by discounting future expenses to present dollar 
values. 

To illustrate the evaluation of P^ and consider a simple economic 
situation in which the only significant costs are fuel and system equip- 
ment costs. The fuel cost is assumed to escalate at a constant annual 
rate, and the owner pays cash for the system. Here, P-j accounts for fuel 
escalation and the discounting of future payments. The factor P^ accounts 
for investment related expenses which in this case, consist only of the 
investment which is already expressed in current dollars. The factors P^ 
and ?2 are then 

P 1 = PWF(N, e, d) (2) 


where N = Period of economic analysis (yrs) 
e = Escalation rate of fuel price 
d = Annual discount rate 

*D1scounting refers to the fact that an expense that is anticipated to be 
$1000 in 10 years is equivalent to an investment today of $463 at a discount 
rate of 8%. 
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The function PWF(N, e, d) Is the present worth factor that accounts for 
inflating payments In discounted money. 


PWF(N, e, d) = — 3 — [. / 1+e \N 

d - e ' \ 1 + d J 


( 3 ) 


When multiplied by a first period cost (which is Inflated at a rate, e, and 
discounted at a rate, d, over N years), the resulting value is the present 
worth life cycle cost. 


In the more complex analysis the expenditures incurred by the additional 
capital investment cause P^ and P 2 t0 take the following form: 

P 1 ■ (1 - Ct) PWF(N, e, d) (4) 

P 2 = P 21 + P 22 " P 23 + P 24 + p 25 " P 26 ' P 27 
where P 21 3 Factor representing the down payment 

P 22 3 Pac t° r representing the life cycle cost 
of the mortgage principal and interest 

P 22 * Factor representing income tax deductions 
for interest payment 

P 24 3 factor representing miscellaneous costs 
(maintenance, insurance, etc) 

p oc = Factor representing net property tax costs 

CD 

P 26 3 f ac ^ or representing straight line depreciation 
tax deduction for commercial installations 

P 2 y ■ Factor representing salvage (commercial installation) 
or resale value (residential installation). 
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2i through P27 are defined as follows: 


The factors P 


P 21 

= D 



P 22 

= (1 - 

D) PWF (N, 0, d)/PWF (N, 0, 

. i) 

P 23 

= (1 - 

D) t jpWF (N, i, d) [i - 

1/PWF (N, 0, 



+ PWF (N, 0, d)/PWF 

(N 0, i)J, 

<cJ- 

CSJ 

= (1 - 

CT) MPWF (N, g, d) 


P 25 

= t (1 

- t) VPWF (N, g, d) 


P 26 

= (Ct /N) PWF (N, 0, d) 


P 27 

= G/(l 

+ d) N 



(6) 

(7) 

1 )] ( 8 ) 

(9) 

(10) 

(11) 

( 12 ) 


where D = Ratio of down payment to the initial investment 

N = Period of analysis (Note that the period of analysis, 
the term of the loan, the depreciation lifetime, and 
the years over which the depreciation deductions con- 
tribute to the analysis are arbitrarily set equal in 
this report). 

d = Discount rate (after tax return on the best 
alternative investment) 

i = Annual mortgage interest rate 

t = Effective income tax rate 

C = Commercial or non-commercial flag (1 or 0 
respectively) 
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M = Ratio of first year miscellaneous costs to 
initial Investment 

g = General inflation rate 

t * Property tax rate based on assessed value 

V = Ratio of assessed value in first year to Initial 
investment 

G » Ratio of salvage or resale value to Initial 
investment 

For a given location, heating load, and economic situation, it is possible 

to optimize the system design variables to yield the maximum life cycle 
savings. The main solar energy system design variable is the collector 
area. The effect of collector area on the life cycle savings is illustrated 
In Figure 4-1 for the four sets of economic conditions. Curve A corresponds 
to an economic scenario in which solar energy cannot compete with the conven- 
tional system. Curve B exhibits a non-zero optimum area, but the conventional 
system is still the most economical. Curve C corresponds to the critical con- 
dition where solar energy can just compete with the conventional system. 

Curve D corresponds to an economic scenario in which the solar energy system 
is the most economical. 

Each curve of Figure 4-1 begins with a negative savings for zero collector 
area. The magnitude of this loss Is C £ , and reflects the presence of solar 
energy system fixed costs in the absence of any fuel savings. As the col- 
lector area increases Curves B, C, and D show increased savings until reach- 
ing a maximum at some optimum collector area. As the collector area is further 
Increased, the fuel savings continue to increase, but the excessive system 

cost forces the life cycle savings of the system to decrease. These collec- 
tor areas at each of the five analysis sites listed In this report have been 
optimized by the f-Chart program analysis technique for the long-term average 
weather conditions and the economic conditions at that site. 
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COLLECTOR AREA 


Figure 4-1 


Life Cycle Savinqs versus Collector Area 
for Four Sets of Economic Conditions 
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4.2 Federal Tax Credits for Solar Energy Systems 


The Federal Government has provided tax incentives that are applicable to 
solar energy systems after 1979. This credit Is 40 percent of the first 
$10,000 spent on solar equipment, or a maximum credit of $4,000. The 
credit Is applied In this analysis by reducing both the collector area 
dependent cost and the cost Independent of the collector area, or con- 
stant solar cost, by an effective credit factor based on the total cost 
of the system. 

As an example of the tax credit computation, assume the collector area 
dependent cost Is $30/Ft based on 100 Ft and the constant solar cost 
Is $900 for a total price of $3900. The effective credit factor Is 0.4 
since the system cost Is less than $10,000. 

Therefore the adjusted costs used as f-Chart Inputs are: 

Collector area dependent cost 

C^'VWx (1 - 0.4) = $18. 00/Ft 2 

Constant solar cost 

C E , = $900 x Cl - 0.4) = $5)40 

If the system cost had exceeded $10,000 the effective credit factor would 
have been the ratio of the maximum credit C$4,000) to the total system cost. 

The f-Chart economic analysis Is modified by using these adjusted costs to 
reflect tax credit effects. Including tax credit In area optimization Is 
an iterative process since the credit is affected by the system size and 
vice versa. Optimal collector area is modified In this analysis, as are 
the f-Chart economic parameters, by use of the tax credit. Items 23 and 24 
in Table 5.1-2 reflect the solar costs after application of tax credits in 
terms of collector area dependent cost and constant cost. Initial system 
costs before and after tax credit inclusion are shown in Table 5.2-1 for 
each site based on optimal collector area. 
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5. RESULTS OF ANALYSIS 


5.1 Technical Results 


For each of the five analysis sites an optimal solar system based on the 
configuration of the actual installation is determined by using the f-Chart 
design procedure. The environmental parameters and the loads used in this 
procedure for each of the five sites are shown in Table 5.1-1. In applying 
the design procedure a process that iterates on the collector area is used. 
Figures 5.1-1 (a) - (e) show the results of that design procedure in terms 
of the expected solar fraction versus the collector area for each site. 

The expected solar fraction is the ratio of the expected solar energy 
used toward satisfying the load to the total load. The graphs in Figures 
5.1-1 (a) - (e) show that as the collector areas increases, the expected 
solar fraction increases. However, the economically optimal collector 
area was selected to maximize the economic benefits of the solar energy 
system, not the expected solar fraction. The optimal collector area is 
shown by the dotted line for each site. Increasing the collector area 
beyond the optimal value forces a diminishing return on the investment for 
the system. The expected solar fraction for the optimal collector area is 
shown in the last column in Table 5.1-1. 

The resulting thermal performance, once the optimal size system is selected, 
is shown in the graphs of Figures 5.1-2 (a) - (e) for each analysis site. 

The incident solar energy is derived from long-term average insolation at 
the site. The total load is computed based on design parameters of the 
actual system as installed, modified by environmental conditions at each 

site. The load calculations are detailed in Appendix A. The useful solar 
energy is the product of the system solar fraction and the total load. It 
shows on a month by month basis the portion of the total load that is ex- 
pected to be supplied by solar energy. The shaded portion between the 
total load curve and the curve of useful solar energy must be supplied by 
conventional energy. 
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As shown in Figures 5.1-1 (a) - (e), the optimal collector areas vary 
from a low of 156 square feet in Fort Worth, Texas to a high of 273 square 

feet at Albuquerque, New Mexico and Washington, D.C. Albuquerque also 
has the highest solar fraction (76%) at the optimal area due to the high 
annual dally insolation and lowest hot water load. (See Table 5.1.1.) 

The technical parameters that describe the solar energy system are listed 

in Table 5.1-2 as Items 1 through 21 and Items 47 and 48 and described in 
detail in Appendix A. Their values are listed by site in Table 5.1-3. The 
remaining technical parameters are assigned values which are constant for 
all sites. 

The economic parameters for the solar energy system are listed in Table 
5.1-2 as Items 22 through 46, and are also described in Appendix A with 
the source for the assigned value designated. 

The following Items are a function of the analysis site. 

• Collector area 

• Collector slope 

• Azimuth angle 

• Effective building UA (applicable to space heating systems) 

• Water main temperature 

• Present cost of solar backup fuel 

• Present cost of conventional fuel 

These are listed by site in Table 5.1-3. 
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The lowest solar fraction (27%) at Akron, Ohio is attributed to the 
high heating and hot water loads and the lowest annual average daily 
insolation of all the analysis sites. In each case the solar fraction 

is proportional to the amount of solar energy available. The optimal 
collector area depends on both the solar energy available and the cost 
of conventional energy. 

Figure 5.1-2 (b) shows that this system can supply almost the entire 

heating and hot water energy requirement from April through November 
for the Albuquerque site. It should also be noted that supplemental 
auxiliary energy is required throughout the year at the Akron, Fort 
Worth, Madison, and Washington, D.C. sites. 
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SUMMARY TABLE 


TABLE 5.1-1 


SOLAR SYSTEM LOAD FACTORS AND ENVIRONMENTAL PARAMETERS 


SITE 

TOTAL ANNUAL LOAD 
(MILLION BTU) 

ENVIRONMENTAL PARAMETERS 

LONG TERM 

EXPECTED 

SOLAR 

FRACTION* 

HEATING 

HOT 

WATER 

INSOLATION 

btu/ft 2 -day 

HEATING 
DEGREE DAYS 

SUPPLY WATER 
TEMP (°F) 

AKRON 

35.44 

26.82 

1086 

6224 

63 

26.6 

ALBUQUERQUE 

26.99 

22.88 

1828 

4292 

73 

75.8 

FORT WORTH 

14.76 

26.23 

1475 

2382 

65 

49.3 

MADISON 

1 

40.30 

30.54 

1 

1191 

7730 

54 

29.9 

WASHINGTON 

30.24 

28.33 

1208 

5010 

60 

42.9 


*For optimal collector area 
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Figure 5.1-1 (b) Solar Fraction vs Collector Area for Albuquerque, New Mexico 
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Figure 5.1-1 (c) Solar Fraction vs Collector Area for Fort Worth, Texas 
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Figure 5.1-1 (d) Solar Fraction vs Collector Area for Madison, Wisconsin 
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Figure 5.1-1 (e) Solar Fraction vs Collector Area for Washington, D. C. 
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Figure 5. 1-2 (a) Thermal Performance of Solar Energy System with Optimized Collector Area for Akron, Ohio 
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Figure 5 . 12 (b) Thermal Performance of Solar Energy System with Optimized Collector Area 
for Albuquerque, New Mexico 
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Figure 5. 1-2 (c) Thermal Performance of Solar Energy System with Optimized Collector Area 
for Fort Worth, Texas 
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Figure 5. 1-2 (d) Thermal Performance of Solar Energy System with Optimized Collector Area 
for Madison , Wisconsin 


32 




MILLION BTU 


20 


WASHINGTON, D. C. 

OPTIMAL COLLECTOR AREA = 273 FT? 


18 


16 


14 


12 


10 


8 


6 


4 


2 


0 



JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC 


MONTH 


Figure 5. 1-2 (e) Thermal Performance of Solar Energy System with Optimized Collector Area 
for Washington, D. C. 
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TABLE 5.1-2 


f-CHART INPUT VARIABLES 


ITEMS VARIABLE DESCRIPTION VALUE 

1 AIR SH+WH = 1, LIQ SH+WH = 2, AIR OR IQ WH ONLY = 3 . . . 1 

2 IF 1, WHAT IS (FLOW RATE/COL. AREA)(SPEC. HEAT)? 2.09 

3 IF 2, WHAT IS (EPSILON)(CMIN)/(UA)? N/A 

*4 COLLECTOR AREA 

5 FRPRIME-TAU-ALPHA* PRODUCT* (NORMAL INCIDENCE) ’. . *. !*.’.! 0.51 

6 FRPRIM-UL PRODUCT 0.65 

7 INCIDENT ANGLE MODIFIER (ZERO IF NOT AVAIL.) 0 

8 NUMBER OF TRANSPARENT COVERS 2 

*9 COLLECTOR SLOPE 

* 10 AZIMUTH ANGLE (E.G. SOUTH = 0, WEST =90) 

11 STORAGE CAPACITY 11.11 

* 12 EFFECTIVE BUILDING UA 

13 CONSTANT DAILY BLDG. HEAT GENERATION 

14 HOT WATER USAGE 138.0 

15 WATER SET TEMP. (TO VARY BY MONT, INPUT NEG.#) 127.0 

16 WATER MAIN TEMP (TO VERY BY MONTH, INPUT NEG. #) 

17 CITY CALL NUMBER 2 

18 THERMAL PRINT OUT BY MONTH = 1 , BY YEAR = 2 1 

19 ECONOMIC ANALYSIS ? YES = 1 , NO = 2 1 

20 USE OPTMZD. COLLECTOR AREA = 1 , SPECFD. AREA =2 

21 SOLAR SYSTEM THERMAL PERFORMANCE DEGRADATION 

22 PERIOD OF THE ECONOMIC ANALYSIS 20 

23 COLLECTOR AREA DEPENDENT SYSTEM COSTS Note 1 

24 CONSTANT SOLAR COSTS Note 1 

25 DOWN PAYMENT (% OF ORIGINAL INVESTMENT) 20 

26 ANNUAL INTEREST RATE ON MORTGAGE 13.5 

27 TERM OF MORTGAGE 20 

28 ANNUAL NOMINAL (MARKET) DISCOUNT RATE 8.5 

29 EXTRA INSUR./MAINT. IN YEAR 1 (% OF ORIG. INV.) 0.5 

30 ANNUAL % INCREASE IN ABOVE EXPENSE 10.0 

* 31 PRESENT COST OF SOLAR BACKUP FUEL (BF) 

32 BF RISE: %/YR = 1, SEQUENCE OF VALUES = 2 1 

33 IF 1, WHAT IS THE ANNUAL RATE OF BF RISE 12.5 

* 34 PRESENT COST OF CONVENTIONAL FUEL (CF ) Note 2 

35 CF RISE: %/YR = 1, SEQUENCE OF VALUES = 2 1 

36 IF 1, WHAT IS THE ANNUAL RATE OF CF RISE 12.5 

37 ECONOMIC PRINT OUT BY YEAR = 1, CUMULATIVE = 2 1 

38 EFFECTIVE FEDERAL - STATE INCOME TAX RATE 30 

39 TRUE PROP. TAX RATE PER $ OF ORIGINAL INVEST 0 


UNITS 


BTU/H-°F-FT 2 
(TABLE 5.1-3) 
BTU/H-°F-FT 2 


(TABLE 5.1.3) 
(TABLE 5.1-3) 
BTU/°F-FT Z 
(TABLE 5.1-3) 
(TABLE 5.1-3) 
GAL/ DAY 
°F 

(TABLE 5.1-3) 


(TABLE 5.1-3) 

YEARS 

$/FT 2 

$ 

% 

% 

YEARS 

% 

% 

% 

(TABLE 5.1-3) 

% 


% 

% 

% 
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TABLE 5.1-2 


f-CHART INPUT VARIABLES (Continued) 


ITEMS VARIABLE DESCRIPTION VALUE UNITS 


40 ANNUAL % INCREASE IN PROPERTY TAX RATE N/A 

41 CAL. RT. OF RETURN ON SOLAR INVTMT? YES = 1, NO = 2 

42 RESALE VALUE (% OF ORIGINAL INVESTMENT) 0 % 

43 INCOME PRODUCING BUILDING? YES = 1, NO = 2 2 

44 DPRC. : STR.LN=1 ,DC.BAL.=2,SM-YR-DGT=3,N0NE=4 2 

45 IF 2, WHAT % OF STR.LN DPRC. RT IS DESIRED? 150 % 

46 USEFUL LIFE FOR DEPREC. PURPOSES 20 YEARS 

47 ECONOMIC COEFFICIENT OF PERFORMANCE OF BACKUP HEATING SYSTEM . TABLE 5.1-3 

48 ECONOMIC COEFFICIENT OF PERFORMANCE OF BACKUP WATER HEATER . . 1 


NOTE: 1. The values of Collector Area Dependent System Costs and Constant Solar 
Costs depend on system size (because of the Federal Tax Credit). These 
costs are listed in Table 5.2-1. The Area Dependent Cost listed in 
Table 5.2-1 must be divided by the optimal area to obtain the value for 
Collector Area Dependent System Costs. 

NOTE: 2. Since the backup for the solar system Is assumed to be the same type 
of system as would conventionally be used without a solar system, 
backup fuel costs and conventional costs per million Btu are equal. 
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TABLE 5.1-3 


SOLAR SYSTEM TECHNICAL PARAMETERS FOR F-CHART PROGRAM 


VARIABLE DESCRIPTION 

UNITS 

LOCATION 

AKRON 

ALBUQUERQUE 

FORT WORTH 

MADISON 

WASHINGTON 

COLLECTOR AREA- OPTIMAL 

FT 2 

195 

273 

156 

195 

273 

COLLECTOR SLOPE 

DEGREES 

45 

45 

43 

53 

49 

AZIMUTH ANGLE 

DEGREES 

0 

0 

0 

0 

0 

EFFECTIVE BLDG UA 

BTU/°F-DAY 

7200 

8081 

8901 

6515 

7745 

CONSTANT DAILY BLDG HEAT GENERATION 

BTU/DAY 

33840 

33840 

33840 

33840 

33840 

SUPPLY WATER TEMPERATURE 

°F 

SEE TABLE l 

C-l FOR MONTHLY 

VALUES 



SYSTEM THERMAL PERF. DEGRADATION 

%/YR 

0 

0 

0 

o 

0 

PRESENT COST OF SOLAR BACKUP FUEL* 

$/MMBTU 

17.50 

20.39 

13.01 

12.21 

19.78 


$/KWH 

0.060 

0.070 

0.044 

0.042 

0.068 

COEFF. OF PERF. FOR HEAT PUMP 


2.26 

2.39 

2.46 

2.25 

2.39 


* An effective rate is computed for each location based on 1000 kWh usage. This effective rate includes all 
charges specified in the rate schedules in Appendix D. 




5.2 Economic Results 


An essential factor in maximizing the life cycle savings of a solar 
energy system, or conversely, of minimizing life cycle costs Is the 
economic optimization of the collector area based on equipment and 
fuel (conventional energy) costs and the capability of the solar sys- 
tem to replace significant quantities of conventional energy with 
solar energy. The replacement capability Is directly dependent on 
the environmental conditions at the Installation site, i.e. available 
solar energy. 

The graphs of Figures 5.2-1 (a) - (e) show the relationship of the factors 
comprising life cycle costs - equipment costs and fuel costs - as a func- 
tion of collector area. Both costs are presented in terms of present 
value, i.e. baselined to today's dollars. It can be readily seen that 
as collector area increases, solar equipment costs increase proportion- 
ately. Also, as collector area increases the fuel costs decrease, 
although not as a straight line function. At some given collector area, 
the sum of these two costs is a minimum, as shown by the life cycle cost 
(LCC) curve. This minimum defines the optimal collector area for the 
given installation site. 

The solar equipment costs discussed in the preceding paragraphs include 
the principal and interest paid on a 13.5 percent, 20 year mortgage, the 
income tax deduction for interest for an owner in the 30 percent bracket 
and the insurance and maintenance costs estimated at 0.5 percent of the 
initial costs. The fuel cost is that which is required by the conven- 
tional backup system and includes the effects of the f-Chart solar system 
model . 

The life cycle costs are not to be confused with life cycle savings. 

Life cycle savings is the difference between the life cycle costs of 
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fuel for a conventional system and the life cycle cost of owning, 
operating and maintaining a solar energy system. 

The life cycle cost (LCC) for each analysis site is clearly defined in 
Figures 5.2.1 (a) - (e). The lowest LCC was found in Fort Worth, 

Texas and the highest LCC was found in Akron, Ohio. Factors which 
have the greatest influence on the LCC are collector area, backup, 
fuel cost and solar fraction. The best combination of these factors 
provide the lowest LCC. 

A summary of the costs and savings for the conventional system and the 
solar energy system is shown in Table 5.2-1 in terms of today's dollars 
expended over the analysis period. It should be recalled that the equip- 
ment costs shown do not include the cost of the conventional system since 
this system must be provided with or without the solar energy system. 

The equipment costs include only the additional hardware that must be 
provided for the solar energy system. This includes the following: 

• Collectors and mounting hardware 

• Piping and duct work (including valves and dampers) 

• Heat exchanger(s) 

• Storage unit(s) 

• Control system 

The best estimates of equipment costs for solar energy systems indicate 
that costs fall into two categories; (1) costs dependent on collector area 
and, (2) costs independent of collector area, or constant costs. This is 
the case, especially for residential systems, because regardless of the 
exact collector area used, certain items of equipment must be provided and 
the costs of hardware and labor for installation seem to be relatively 
constant. However, the cost of collectors, and certain incremental costs, 
are dependent on the size of the collectors used. These costs are shown 
in Table 5.2-1 for each of the five analysis sites and the total cost for 
the system is the sum of the constant cost and the area dependent cost 
multiplied by the collector area. 
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The Initial cost of the system In this analysis should be adjusted 
for the federal tax credit (and any other tax credit allowed by the 
state or local governments) by the methods discussed In Section 4.2. 

These adjusted costs are shown In parentheses under "Initial Cost of 
System" In Table 5.2-1 and are used In computing the "Present Worth 
of Total Solar Costs." 

Some conventional energy must be expended with or without the solar 
energy system because. In most cases, the solar energy system will 
replace only a portion of the total energy required to support the 
load. Savings are possible with the solar system only because the 
total costs with the solar system are less than the costs of conven- 
tional energy. Consequently, the fuel costs over the analysis period 
(20 years) are shown In Table 5.2-1 with and without the solar system. 

It Is assumed In this analysis that the solar system would be financed 
through a 20 year loan at an Interest rate of 13.5 percent. Property taxes 
are assumed to be zero, but this may not be universally true. Insurance 
on the value of the solar energy system and maintenance costs are assumed 
to be 0.5 percent per year of the initial costs. Since Interest paid on 
a loan Is tax deductible for federal taxes, and in most cases for state 
taxes, at different rates according to the income tax bracket of the 
borrower, a 30 percent combined federal -state tax bracket was assumed. 

The value of all these costs based on the assumptions of this analysis 
is shown as the "Present Worth of Other Solar Costs" in Table 5.2-1. 
Combined with the costs for fuel with the solar system, the value is 
the "Present Worth of Total Solar Costs." 

Since only incremental equipment and associated costs are Included In the 
analysis, the present worth of total costs for the conventional system 
without solar are simply the cost of fuel without solar. Then the "Present 
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Worth of Cumulative Savings" is the difference between the "Present 
Worth of Total Costs Without Solar" and the "Present Worth of the Total 
Costs With Solar". These values for each of the five analysis sites 
are listed in Table 5.2-1. 

Finally, two economic performance parameters called "Year of Positive 
Savings" and the "Year of Payback" are shown in Table 5.2-1. As previ- 
ously discussed the year of positive savings is the year after purchase 
in which the solar system first becomes profitable, i.e., the annual 
conventional fuel bill without solar exceeds sum of the annual fuel bill 
with solar and the annual costs for the solar system. The year of payback 
is the year after purchase when the compounded net savings equals the 
initial cost for the solar system. Savings are compounded at the discount 
rate throughout the analysis period. The factors that determine years until 
positive savings are shown in Figures 5.2.2 (a) - (e) for each analysis site. 
The factors that determine the years until payback are shown in Figures 5.2-3 
(a) - (e) for each analysis site. The year corresponding to the intersection 
of the "Mortgage Principle Remaining" curve and the "Compounded Solar Savings" 
curve is the year that the savings are sufficient to pay off the mortgage 
balance. 

As shown in Table 5.2-1, a solar energy system of the type installed at the 
Akron, Ohio site is not economically feasible for the regions of Akron, 

Fort Worth, Madison, and Washington, D.C. Figures 5.2-2 (a) - (e) graphically 
illustrate that a positive savings occurs on each site beginning in five 
years at Albuquerque, N.M., eight years at Washington, D.C., eleven years at 
Fort Worth, Texas, twelve years at Akron, Ohio, and thirteen years at Madison, 
Wisconsin. Conventional energy cost actually dictates the occurrence of 
"positive savings" for each site. 

The "Year of Payback" is graphically illustrated in Figures 5.2-3 (a) - (e) 
and show that only the Albuquerque, New Mexico site had a reasonable payback 
period of less than twenty years (actually seventeen years). All other payback 
periods exceeded the reasonable twenty year and basically the conventional 
energy cost dictated the occurrence of payback. 


40 



PRESENT WORTH IN THOUSANDS OF DOLLARS 


AKRON, OHIO . 

OPTIMAL COLLECTOR AREA = 195 FT 2 



COLLECTOR AREA (FT 2 ) 


Figure 5.2-1 (a) Optimization of Collector Area For Akron, Ohio 









PRESENT WORTH IN THOUSANDS OF DOLLARS 


ALBUQUERQUE, NEW MEXICO , 

OPTIMAL COLLECTOR AREA - 273 FT 2 
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Figure 5J-1 (b) Optimization of Collector Area for Albuquerque, Nm Mexico 
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Figure 5.2-1(c) Optimization of Collector Area for Fort Worth , Texas 






PRESENT WORTH IN THOUSANDS OF DOLLARS 


MADISON, WISCONSIN 

OPTIMAL COLLECTOR AREA = 195 FT2 



COLLECTOR AREA ( FT2 ) 


Figure 5 2-1 (d) Optimization of Collector Area for Madison, Wisconsin 



PRESENT WORTH IN THOUSANDS OF DOLLARS 


WASHINGTON, D. C 

OPTIMAL COLLECTOR AREA = 273 FT2 
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Figure 5.2-1 (e) Optimization of Collector Area for Washington, D. C. 




SUMMARY TABLE 


TABLE 5.2-1 

COSTS AND SAVINGS OVER 20 YEAR ANALYSIS PERIOD IN DOLLARS (1980) 



INITIAL COST OF SYSTEM 1 

PRESENT WORTH 
OF FUEL COSTS 

PRESENT 

WORTH 

OF 

OTHER 

PRESENT 

WORTH 

OF 

TOTAL 

PRESENT 
WORTH 
OF TOTAL 

PRESENT 

WORTH 

OF 

YEAR OF 


SITE 

CONSTANT 

ARFfl 

DEPENDENT 

TOTAL 

gg 




COSTS W/0 
SOLAR 

CUMULATIVE 

SAVINGS 

POSITIVE 

SAVINGS 

YEAR OF 
PAYBACK 

AKRON 

8519 

(5623) 

, r 

3231 

(2133) 

11750 

(7756) 

13657 

19754 

9033 

22690 

19764 

-2926 

12 

>20 

ALBUQUERQUE 

8519 

(5904) 

4524 

(3132) 

13043 

(9038) 

4171 

18516 

10539 

14710 

18516 

3805 

5 

17 

; FORT WORTH 

8519 

(5452) 

2595 

(1654) 

11114 

(7106) 

5181 

11141 

8280 

13461 

11141 

i 

-2320 

11 

>20 

MADISON 

8519 

(5623) 

3231 

(2133) 

11750 

(7756) 

10415 

15719 

9046 

19461 

15719 

-3741 

13 

>20 

1 WASHINGTON 

t 

8519 

(5904) 

4524 

(3134 

13043 

(9038) 

11292 

1 

j 

21539 
.1 

10535 

21827 

21539 

1 

-288 

8 

i 

>20 


NOTE: 


1. Values in parentheses are sdjusted for the Federal tax credit by the method detailed in Section 4,2. 











ANNUAL COST IN DOLLARS 



YEARS 


Figure 5.2-2 (a) Annual Expenses for Solar and Conventional System at Akron, Ohio 
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ALBUQUERQUE, NEW MEXICO 
OPTIMAL COLLECTOR AREA = 273 Fr 



YEARS 


Figure 5 2-2 (b) Annual Expenses for Solar and Conventional System at Albuquerque, New Mexico 
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ANNUAL COST IN DOLLARS 



YEARS 

Figure 5.2-2(c) Annual Expenses for Solar System and Conventional System at Fort Worth, Texas 
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ANNUAL COST IN DOLLARS 


MADISON, WISCONSIN 

OPTIMAL COLLECTOR AREA = 195 FT2 
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Figure 5J-2(d) Annual Expenses for Solar System and Conventional System at Madison, Wisconsin 
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ANNUAL COST IN DOLLARS 


WASHINGTON, DC . 

OPTIMAL COLLECTOR AREA = 273 Fr 



YEARS 

Figure 52-2 fe) Annual Expenses for Solar and Conventional System at Washington, D.C. 
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Figure 5.2-3 (a) Payback for Solar Energy System for Akron, Ohio 
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ALBUQUERQUE, NEW MEXICO 
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Figure 5.2-3 (b) Payback for Solar Energy System for Albuquerque, New Mexico 
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Figure 5.2 3 (c) Payback for Solar Energy System for Fort Worth, Texas 
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DOLLARS IN THOUSANDS 
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Figure 5.2-3 (d) Payback for Solar Energy System for Madison, Wisconsin 
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YEARS 

Figure 5.2-3 (e) Payback for Solar Energy System for Washington , D.C. 
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6. ECONOMIC UNCERTAINTY ANALYSIS 


The economic evaluation methods presented In this report are based on 
the assumption that reliable values for economic variables can be as- 
signed. However, there Is an inherent uncertainty In predicting future 
expenses and benefits which Is magnified by international economic unsta- 
bility. As a consequence, the results of both the life cycle cost analysis 
and the optimization procedures must be accepted with discretion and the 
effect of uncertainties must be evaluated. 


For a given set of conditions, the change In the present worth of life 
cycle cumulative savings (Table 5.2-1), ALCCS, resulting from a change In 
a particular variable. Ax., can be approximated by the following: 

J 


. aLCCS 

4LCCS ' ST. 4x j 

J 


(13) 


The expression for 9LCCS/3X. can be obtained by direct differentiation of 

J 

the life cycle savings equation. The life cycle cost model of Equations 
(1), (4) and (6)-(12) will be used for this analysis. The derivatives of 
these equations for each variable are given in Appendix B. To Illustrate 
the use of these relationships. Uncertainty Analysis Tables 6-1 through 
6-5 were made up for each analysis site. The tables give the change in 
solar system life cycle cumulative savings, ALCCS, caused by a 10 percent 
relative increase In each of the variables. 


Table 6-1 shows, for example, that a 10 percent increase in the discount 
rate from 8.5 to 9.4 percent yields a decrease in the value of of 
approximately 2.43 giving a modified value of = 24.14. The value of 
?2 decreases by 0.065 giving a modified value of P 2 = 1.099. The value of 
LCCS decreases by approximately $21 or a relative change of 1 percent 
in the baseline value of $2,926. By comparing the magnitude of ALCCS for 
each variable, the relative sensitivity of the savings to a change In the 
variable can be assessed. From the table. It is evident that the savings 
are affected most by a change in fuel inflation, and least by a change in 
down payment. The complex relationship of the variables to each other makes 
an intuitive approach unreliable and necessitates analysis of this type. 
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The information of Tables 6-1 through 6-5 can also be used to estimate the 
uncertainty in life cycle cumulative savings due to uncertainty in different 
variables. If all the economic parameters are subject to variation 
a reasonable estimate of savings uncertainty can be obtained by the 
following: 


ALCCS. 


= 

£ ( 

3LCCS 

)■] 


^ \ 

3x. 


L j = i 

J 



As an example, assume uncertainties of +10 percent in all seventeen of the 
variables listed in Table 6-1. The probable uncertainty estimate, using 
the data from the Table is: 



ALCCS prob = $2935 


Cumulative Savings = $3805 


It should be noted that the ALCCS value is less than the cumulative savings 
(See Table 5.2-1) for the Albuquerque, New Mexico site. For a reasonable 
and favorable change in all the economic variables listed in Table 6-1, there 
is little possibility of a savings with this system. The results for other 
sites are as follows: 


Akron, Ohio 

ALCCS prob = $1324 
Cummulative Savings = -$2926 

Fort Worth, Texas 

ALCCS prob = $1349 
Cumulative Savings = -$2320 

Madison, Wisconsin 

ALCCS prob = $1247 
Cumulative Savings = -$3741 
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ALCCS prob = $2037 
Cumulative Savings * -$288 
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TABLE 6-1 

UNCERTAINTY ANALYSIS FOR AKRON, OHIO 


2 

Optimized Collector Area = 195 FT 


COST PARAMETER (x.) 

NOMINAL 

VALUES 

NOMINAL 

VALUE 

DELTA 

3P1 

3X j 

3P2 

ax j 

3LCCS 

3x j 

ALCCS 

AREA DEPENDENT COST (c a ) 

T03 

1.0930 

0.0 

0.0 

1227 

=248 

AREA INDEPENDENT COST (C £ ) 

5619.00 

516.9000 

0.0 

0.0 

-1 

-654 

ELECTRICAL ENERGY COST (C fE ) 

17.50 

1.7500 

0.0 

0.0 

301 

1 526 

DOWN PAYMENT/ I NIT INV. (D) 

0.200 

0.0200 

0.0 

-0.074 

571 

11 

FIRST YR. MISC COST/INIT INV (M) 

0.005 

0.0005 

0.0 

21.066 

-163267 

-82 

FIRST YR. ASSESSED VAL/INIT INV (V) 

0.0 

0.0 

0.0 

0.0 

0 

0 

SALVAGE VAL/INIT INV (G) 

0.0 

0.0 

0.0 

mm 

1516 

0 

ANNUAL MKT DISCOUNT RATE (d) 

0.085 

0.0085 

-286.35 


2453 

21 

ANNUAL MKT RATE OF FUEL COST INC. (e) 

0.125 

0.0125 

252.55 

MM 

49965 

625 

ANNUAL INT. RATE ON MORTGAGE (i) 

0.135 

0.0135 

0.0 

4.406 

-34150 

-461 

ANNUAL RATE OF GENERAL INFLATION (g) 

0.100 

0.0100 

0.0 

0.954 

-7392 

-74 

PROPERTY TAX RATE (t) 

0.0 

0.0 

0.0 

0.0 

0 

0 

EFFECTIVE INCOME TAX RATE (t) 

0.300 

0.0300 

0.0 

-0.838 

6492 

195 

ANNUAL HOT WATER LOAD (L £ ) 

26.82 

2.6820 

0.0 

0.0 

124 

332 

ANNUAL HEATING LOAD (L f ) 

35.44 

3.5440 

0.0 

0.0 

55 

194 

ANNUAL SOLAR FRACTION (F) 

0.266 

0.0266 

0.0 

0.0 

19762 

526 

COEFFICIENT OF PERFORMANCE (COP p ) 

2.260 

0.2260 

0.0 

0.0 

-858 

-194 



TABLE 6-2 


UNCERTAINTY ANALYSIS FOR ALBUQUERQUE, NEW MEXICO 


Optimized Collector Area » 273 FT^ 


COST PARAMETER (Xj) 

NOMINAL 

VALUES 

NOMINAL 

VALUE 

DELTA 

3P1 

3X J 

3P2 

3X J 

3LCCS 

3X . 

j 

ALCCS 

AREA DEPENDENT COST (C A ) 

11.49 

1.1490 

0.0 

0.0 

-318 

-365 

AREA INDEPENDENT COST (C p ) 

5906.00 

590.6000 

0.0 

0.0 

-1 

-688 

ELECTRICAL ENERGY COST (C pE ) 

20.39 

2.0390 

0.0 

0.0 

688 

1404 

DOWN PAYMENT/ IN IT INV. (D) 

0.200 

0.0200 

0.0 

-0.074 

666 

13 

FIRST YR. MISC COST/INIT INV (M) 

0.005 

0.0005 

0.0 

21.066 

-190492 

-95 

FIRST YR. ASSESSED VAL/INIT INV (V) 

0.0 

0.0 

0.0 

0.0 

0 

0 

SALVAGE VAL/INIT INV (G) 

0.0 

0.0 

0.0 

-0.196 

1769 

0 

ANNUAL Mia DISCOUNT RATE (d) 

0.085 

0.0085 

-286.35 

-7.626 

-82277 

-699 

ANNUAL MKT RATE OF FUEL COST INC. (e) 

0.125 

0.0125 

252.55 

0.0 

133387 

1667 

ANNUAL INT. RATE ON MORTGAGE (i) 

0.135 

0.0135 

0.0 

4.406 

-39845 

-538 

ANNUAL RATE OF GENERAL INFLATION (g) 

0.100 

0.0100 

0.0 

0.954 

-8625 

-86 

PROPERTY TAX RATE (t) 

0.0 

0.0 

0.0 

0.0 

0 

0 

EFFECTIVE INCOME TAX RATE (t) 

0.300 

0.0300 

0.0 

-0.838 

7575 

227 

ANNUAL HOT WATER LOAD (L £ ) 

22.88 

2.2880 

0.0 

0.0 

411 

940 

ANNUAL HEATING LOAD (L p ) 

26.99 

2.6990 

0.0 

0.0 

172 

464 

ANNUAL SOLAR FRACTION (F) 

0.758 

0.0758 

0.0 

0.0 

18514 

1403 

COEFFICIENT OF PERFORMANCE (COP p ) 

2.390 

0.2390 

0.0 

0.0 ! 

-1940 

-464 



TABLE 6-3 


UNCERTAINTY ANALYSIS FOR FORT WORTH, TEXAS 


2 

Optimized Collector Area = 156 FT 


COST PARAMETER (xj 

J 

NOMINAL 

VALUES 

NOMINAL 

VALUE 

DELTA 

3P1 

3X J 

3P2 

3x o 

3LCCS 

sx j 

ALCCS 

Area dependent cost (c a ) 1 

10.60 r 

1.0600 

0.0 

“O 

~ -T82 


AREA INDEPENDENT COST (C E ) 

5450.00 

545.0000 

0.0 

0.0 

-1 

-635 

ELECTRICAL ENERGY COST (C pE ) 

13.01 

1.3010 

0.0 

0.0 

423 

549 

DOWN PAYMENT/ IN IT INV. (D) 

0.200 

0.0200 

0.0 

-0.074 

523 

10 

FIRST YR. MISC COST/INIT INV (M) 

0.005 

0.0005 

0.0 

21.066 

-149642 

-75 

FIRST YR. ASSESSED VAL/INIT INV (V) 

0.0 

0.0 

0.0 

0.0 

0 

0 

SALVAGE VAL/INIT INV (G) 

0.0 

0.0 

0.0 

-0.196 

1390 

0 

ANNUAL MKT DISCOUNT RATE (d) 

0.085 

0.0085 

-286.35 

-7.626 

-5021 

-43 

ANNUAL MKT RATE OF FUEL COST INC. (e) 

0.125 

0.0125 

252.55 

0.0 

52207 

653 

ANNUAL INT. RATE ON MORTGAGE (i) 

0.135 

0.0135 

0.0 

4.406 

-31301 

-423 

ANNUAL RATE OF GENERAL INFLATION (g) 

0.100 

0.0100 

0.0 

0.954 

-6775 

-68 

PROPERTY TAX RATE (t) 

0.0 

0.0 

0.0 

0.0 

0 

0 

EFFECTIVE INCOME TAX RATE (t) 

0.300 

0.0300 

0.0 

-0.838 

5950 

179 

ANNUAL HOT WATER LOAD (L £ ) 

26.230 

2.6230 

0.0 

0.0 

170 

447 

ANNUAL HEATING LOAD (L p ) 

14.760 

1.4760 

0.0 

0.0 

69 

102 

ANNUAL SOLAR FRACTION (F) 

0.493 

0.0493 

0.0 

0.0 

11141 

549 

COEFFICIENT OF PERFORMANCE (COP p ) 

2.460 

0.2460 

0.0 

0.0 

-416 

-102 




TABLE 6-4 


UNCERTAINTY ANALYSIS FOR MADISON, WISCONSIN 


2 

Optimized Collector Area = 195 FT 


COST PARAMETER (x.) 

J 

NOMINAL 

VALUES 

NOMINAL 

VALUE 

DELTA 

3P1 

3X J 

3P2 

SX J 

3LCCS 

3x i 

ALCCS 

AREA DEPENDENT COST (C ft ) 

T 0793 

1.0930 

0.0 

0.0 

-228 

^48 

AREA INDEPENDENT COST (C £ ) 

5619.00 

561.9000 

0.0 

0.0 

-1 

-654 

ELECTRICAL ENERGY COST (C pp ) 

12.21 

1.2210 

0.0 

0.0 

385 

470 

DOWN PAYMENT/ INIT INV. (D) 

0.200 

0.0200 

0.0 

- 0.074 

571 

11 

FIRST YR. MISC COST/INIT INV (M) 

0.005 

0.0005 

0.0 

21.066 

-163267 

-82 

FIRST YR. ASSESSED VAL/INIT INV (V) 

0.0 

0.0 

0.0 

0.0 

0 

0 

SALVAGE VAL/INIT INV (G) 

0.0 

0.0 

0.0 

- 0.196 

1516 

0 

ANNUAL MKT DISCOUNT RATE (d) 

0.085 

0.0085 

- 286.35 

- 7.626 

8454 

72 

ANNUAL MKT RATE OF FUEL COST INC. (e) 

0.125 

0.0125 

252.55 

0.0 

44672 

558 

ANNUAL INT. RATE ON MORTGAGE (i) 

0.135 

0.0135 

0.0 

4.406 

-34150 

-461 

ANNUAL RATE OF GENERAL INFLATION (g) 

0.100 

0.0100 

0.0 

0.954 

-7392 

-74 

PROPERTY TAX RATE (t) 

0.0 

0.0 

0.0 

0.0 

0 

0 

EFFECTIVE INCOME TAX RATE (t) 

0. 300 

0.0300 

0.0 

- 0.838 

6492 

195 

ANNUAL HOT WATER LOAD (L £ ) 

30.54 

3.0540 

0.0 

0.0 

97 

296 

ANNUAL HEATING LOAD (L p ) 

40.30 

4.0300 

0.0 

0.0 

43 

174 

ANNUAL SOLAR FRACTION (F) 

0.299 

0.0299 

0.0 

0.0 

15719 

470 

COEFFICIENT OF PERFORMANCE (COP p ) 

2.250 

0.2250 

0.0 

0.0 

- 772 

-174 




TABLE 6-5 


UNCERTAINTY ANALYSIS FOR WASHINGTON, DC 


Optimized Collector Area 88 273 FT^ 


COST PARAMETER (x.) 

J 

NOMINAL 

VALUES 

NOMINAL 

VALUE 

DELTA 

3P1 

dX . 
J 

3P2 

3X > 
J 

3LCCS 

ALCCS 

AREA DEPENDENT COST (C A ) 

11.49 

1.1490 

0.0 

0.0 

-318 

-365 

AREA INDEPENDENT COST (C £ ) 

5906.00 

590.6000 

0.0 

0.0 

-1 

-688 

ELECTRICAL ENERGY COST (C p£ ) 

19.78 

1.9780 

0.0 

0.0 

46 

924 

DOWN PAYMENT/ I NIT INV. (D) 

0.200 

0.0200 

0.0 

-0.074 

666 

13 

FIRST YR. MISC COST/INIT INV (M) 

0.005 

0.0005 

0.0 

21.066 

-190492 

-95 

FIRST YR. ASSESSED VAL/INIT INV (V) 

0.0 

0.0 

0.0 

0.0 

0 

0 

SALVAGE VAL/INIT INV (G) 

0.0 

0.0 

0.0 

-0.196 

1769 

0 

ANNUAL MKT DISCOUNT RATE (d) 

0.085 

0.0085 

-286.35 

-7.626 

-30620 

-260 

ANNUAL MKT RATE OF FUEL COST INC. (e) 

0.125 

0.0125 

252.55 

0.0 

87827 

1098 

ANNUAL I NT. RATE ON MORTGAGE (i) i 

0.135 

0.0135 

0.0 

4.406 

-39845 

-538 

ANNUAL RATE OF GENERAL INFLATION (g) 

0.100 

0.0100 

0.0 

0.954 

-8625 

-86 

PROPERTY TAX RATE (t) 

0.0 

0.0 

0.0 

0.0 

0 

0 

EFFECTIVE INCOME TAX RATE (t) 

0.300 

0.0300 

0.0 

-0.838 

7575 

227 

ANNUAL HOT WATER LOAD (L £ ) 

28.330 

2.8330 

0.0 

0.0 

225 

639 

ANNUAL HEATING LOAD (L p ) 

30.240 

3.0240 

0.0 

0.0 

94 

285 

ANNUAL SOLAR FRACTION (F) 

0.429 

0.0429 

0.0 

0.0 

21539 

924 1 

COEFFICIENT OF PERFORMANCE (COP p ) 

2.390 

0.2390 

0.0 

0.0 

-1194 

-285 


7. SUMMARY AND CONCLUSIONS 


Solar energy systems of the type installed at the Akron, Ohio site are 
not economically beneficial under the assumed economic conditions at 
Akron, Ohio; Fort Worth, Texas; Madison, WI; and Washington, DC as shown 
in Figure 7-1. Only in Albuquerque, New Mexico, where the average solar 
insolation is 1828 Btu/Ft /day and the conventional energy (electricity) 
cost is high (0.070 $/kWh), is this solar energy system marginally 
profitable. Economic benefits from this solar energy system depend pri- 
marily on two factors: (1) Decreasing the initial investment required; 

(2) The continuing increase in the cost of conventional energy. The 
capability to decrease the cost of the system relative to its present 
level is uncertain. It depends on favorable tax treatment from the 
various levels of government, local through federal, as well as the con- 
tinuing development of the solar energy industry. On the other hand, 
increases in the cost of conventional energy are virtually assured. From 
the economic uncertainty analysis in Section 6, where the conventional 
energy costs are medium to high, the savings with this system are 0.75 
to 1.9 times more sensitive to increases in the solar energy system cost 
than to proportional increases in the conventional energy cost. This 
sensitivity demonstrates the benefit of reducing the investment for the 
solar energy system. 

The analysis and results given in this report can be used to guide a potential 
solar energy system buyer in evaluating the purchase of this type of heating/ 
DHW system. To do this the solar insolation in the buyer's geographic area 
must be known. This data is available from several sources, including [9], and 
[10]. The cost of conventional energy must also be known. The local utility 
company can furnish rates from which a comparison cost based on 1000 kWh 
use can be computed in dollars per kWh. These values can then be compared 
with the characteristics of the analysis sites given in Section 3.1. The 
results for that analysis site can be ascertained from Section 5.1 and 5.2. 

The primary economic parameters such as solar system cost, mortgage rates, 
inflation rates, discount rates, etc., are generally known by the buyer 
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Figure 7-1 Economic Summary Chart for all Analysis Sites 
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for his area. Deviations in these economic parameters from the values 
assumed in developing the results in this report can be evaluated from 
material included in Section 6. The ALCCS values given in Tables 6-1 
through 6-5 were computed based on a 10 percent increase in the economic 
parameter in question. A 10 percent decrease simply means changing the sign 
of the value in the appropriate table. Larger increases or decreases in 
an economic parameter can also be obtained by multiplying the ALCCS value 
by the ratio of the desired increase to the 10 percent increase used 
in the original computation. 

As an example of the discussion above, assume the buyer has determined 
that the characteristics of his locale are similar to Fort Worth, Texas, 
and is considering the results reported for this solar energy system in 
Fort Worth. He notes that the reported loss from Table 5.2-1 Is -$2320; 
however, the conventional energy cost of his locale is $0. 040/kWh, instead 
of the $0. 44/kWh (Table 5.1-3) used in developing the Fort Worth value. 

To modify the value to consider the new rate the change is computed as: 

~ 0 , 4 ~ o To 4 ~ 4 ' 0 ' 4 ' 4 X 100% = 9,1% ( decrease ) 

In Table 6-3 for Fort Worth it can be seen that a 10 percent increase in 
electrical energy cost yields a value for ALCCS of $447. The impact on 
the Life Cycle Cost Savings of a 9.1 percent decrease In fuel cost can 
be computed as follows: 

ALCCS = * $447 = - $407 (decrease) 

Therefore, the new loss is: 

-$2320 - $407 = - $2727 

Consequently the solar energy system moves to an even less competitive 
position because of the lower conventional energy cost. 

The buyer can evaluate the result of a change in any of the economic 
parameters in the same manner. However, he should be aware that the 
parameters are sometimes inter-related and a change in one parameter 
may affect the ALCCS for several parameters. Consequently, the larger 
the changes the less the accuracy. However, approximate results may 
be obtained that prove of value in making a final decision. 
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A-l 



APPENDIX A 


F-Chart Procedure 

Modifications are made to f-Chart to enable the program to be used to 
perform economic analysis of the following: 

1. Systems that use heat pumps and fossil fuel space 
heating systems, as well as electric resistance heat. 

2. Systems that use two different energy sources for 
domestic hot water heating and space heating. 

The problem of analysis of the solar energy system with a conventional 
backup other than electric resistance heat is resolved by Introducing 
Coefficients of Performance (COP's) (Item Nos. 47 and 48) whose values 
are dependent upon the types of backup systems. Typical COP's of heat 
pumps are computed from a heat pump model which uses as inputs the ambient 
and building temperature. Fossil fuel furnace COP's are assumed to be 0.60 
unless different efficiencies, based on manufacturer's or other sources of 
data, are available. 

The problem of analysis with two different energy sources is resolved 
by adjusting the COP's of the space heating system and domestic hot water 
system relative to the cost of electrical energy. This is necessary be- 
cause the structure of f-Chart assumes electric energy to be the source 
for both space heating and domestic hot water. The adjustment factors 
are the adjusted ratios of the rates for the two energy sources used. 

The general expression for this is: 

SH COP' = E1ectrica1 Energy Rate ($/rnmion Btu) , SH cop , 

or ( SH Auxiliary Fuel Rate)($/mill1on Btu) ! or > 

HW COP' j or [ (hWCOP) 

( HW Auxiliary Fuel Rate) 

where the Electrical Energy Rate is the effective rate for 1000 kWh 
and the SH or HW Auxiliary Fuel Rate is the actual cost for fuel 
converted to $/mi 11 ion Btu. Electrical Energy Rate will also be 
used for the value of Items Number 31 and 34 for systems of this 
configuration. 


The value of SH COP' Is Input to the modified f-Chart program. 

This value Is used to compute an adjusted total load. The load, 
in turn. Is used to derive the solar fraction which is input to 
the f-Chart economic analysis subroutine. 

Major considerations of the final report analysis procedure are the 
definitions of the loads that the system supports as it is analyzed 
in different geographic locations, and the sizing of the system to 
handle these loads at the various locations. The method Is outlined 
in the following paragraphs. 

The monthly long-term heating load at the selected analysis sites is 
computed in the f-Chart procedure from the following equation: 

HL LT = UA*HDD lt - HTGEN DAYS 


where 


UA is the modified building energy loss coefficient 

HDD^j is the monthly long-term average heating degree days 

HTGEN is the internally generated heat computed from 
measured data. 

It Is to be noted that UA Is a modified parameter. The modification is 
to compensate for the fact that housing standards differ from location 
to location, i.e. , the construction standards for a Florida house are not 
suitable for the New York environment. The UA factor used Is derived from 
the ASHRAE 90-75 Standard [8] as a function of long term heating degree 
days according to the appropriate U-value. The area. A, is derived from 
the building where the system is Installed. 
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HTGEN Is a factor that accounts for the part of the load which Is 
internally generated. This is assumed to be the heat added which 
brings the building to the desired (comfortable) temperature when 
the outside ambient temperature is 65°F and no auxiliary heat Is 
being added to the building. HTGEN, once derived, Is assumed to 
be constant since it is a function of the life style of the occupants. 
The value of HL L y is the monthly long-term average heat load Input 
to f-Chart. 

Additional technical and economic parameters that are input to f-Chart 
for the final report analysis are listed below with applicable 
comments. 


1. Air SH + WH = 1, Liq SH + WH = 2, Air or Liq WH Only = 3 

Comment: This is a definition of system type. The value 
is 1, if the system uses air collectors and supplies both 
space heat and domestic hot water; 2, if the system uses 
liquid collectors and supplies both space heat and domestic 
hot water; 3, if the system uses either type of collector 
and supplies only domestic hot water. 

2. (Flow rate/ col. area) * (Spec, heat) 

Comment: If the system is an air system, this parameter is 

applicable. It is the air mass flow rate in lb/min divided 
by the gross collector area multiplied by the specific heat 
of air at standard conditions. The value of this parameter 
is computed for the system at the actual installation site. 

This value is then maintained constant as the collector size 
is optimized for all analysis sites.* 

2 

*f-Chart uses an optimized value of 2.15 Btu/Hr-°F-Ft for this parameter. 
In resizing a system, only the collector size is varied. The system is 
not given the benefit of further optimization. 
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3. 


eCmin/UA 


Comment: If the system is a liquid system and uses a liquid 

to air heat exchanger in the space heating loop, this parameter 
is applicable. It is the manufacturer's heat exchanger effec- 
tiveness multiplied by the minimum capacitance rate through 
the heat exchanger and divided by the building energy loss 
coefficient. If the heat exchanger effectiveness is unknown, 
a default value of 0.5 Is specified. The capacitance, Cmin, 
is the- minimum product of mass flow rate and specific heat, 
which usually occurs on the air side. The UA value is the mod- 
ified parameter applicable to the site. Deriving this value 
of UA has been previously discussed. The value of eCmin/UA 
is computed for the system at the actual installation site. 

This value is then maintained constant as the collector size 
is optimized for all analysis sites.* 

4. Collector Area 

Comment: This is the gross collector area which Is optimized 
for all analysis sites. The optimization is extended to the 
actual installation site if an optimum sizing is not apparent 
in the original design. The predicted performance with optimal 
collector sizing is then compared to the predicted performance 
of the actual design and the actual measured performance. 

5. F r (to) 

Comment: The basic value of F R (xa) is derived from the col- 

lector analysis program. This value is more consistent with 
actual operation than the manufacturer's or laboratory single 


*f-Chart uses an optimized value of 2.0 (dimensionless) for this parameter 
In resizing a system only the collector size is varied. 

The system is not given the benefit of further optimization. 
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panel test values. If the system has a heat exchanger 
between collectors and storage, the derived value of 
F r (tcx) was modified by the F R '/F R factor as outlined 
In Section 2.4.4 of EES Report 49-3 (f-Chart Users 
Manual). [4] Note that the values Input to f-Chart are 
assumed to be derived In accordance with ASHRAE specified 
method. 

6. F R U L 

Comment: Same comment as Item 5. 

7. Incidence Angle Modifier 

Comment: In general, the default value of 0 is used. For 

evacuated tube collectors modeled as flat plate collectors 
the collector angle Incidence modifier is obtained from the 
collector manufacturer. 

8. Number of Transparent Covers 

Comment: This is specified according to the characteristics 
of the collector. 

9. Collector Slope 

Comment: Collector Slope Is changed according to the 

latitude of the site and the type of system. When the site 
analyzed is the existing site, the actual slope value is 
used. For other analysis sites the slope Is computed as 
f ol 1 ows : 

t Latitude +10° If space heat and domestic hot water 
• Latitude if domestic hot water only 
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10. Azimuth Angle 

Comment: At sites other than the existing Installation site the 
azimuth angle is 0°. At the existing site the actual azimuth 
angle was used for analysis. However, any resulting performance 
degradation Is noted. 

11. Storage Capacity 

Comment: This parameter Is computed as the product of storage 
mass and specific heat divided by collector area for the exis- 
ting site. The same value of storage capacity Is used for all 
sites. 

12. Effective Building UA 

Comment: The building UA, if not known, is derived from the 
measurement data contained In the Seasonal Report [3]. The 
computed value of UA Is compared for reasonableness with a 
corresponding value of UA derived from ASHRAE Standard 90-75. 

For other analysis sites the value of UA Is derived from 
ASHRAE 90-75 as a function of building type and heating 
degree-days for each site. 

13. Constant Daily Building Heat Generation 

Comment: For residential type buildings, this parameter Is 

derived from the measurement data contained In the Seasonal 
Report [3]. The derived value is held constant for all analysis 
sites. 

14. Hot Water Usage 

Comment: An effective average hot water consumption rate 

that accounts for actual load plus standby losses was 
computed from the following equation: 
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HWCSMPEFF = g- j TOl/oSET j *' TTSET - TMA!NT'*~MS j TMAlN~ + "TSET j 
Number of Days In Month 

15. Water Set Temperature 

Comment: The actual value of this parameter at the existing site 
Is used for all analysis sites. 

16. Water Main Temperature 

Comment: The Inputs for this parameter are a series of monthly 
values. The actual monthly value at the existing site is 
referenced to the average long-term ambient for the month for 
analysis at that site. For analysis at other sites the 
monthly value of TMAIN was established by site measurement 
at a nearby site referenced to the average long-term ambient 
for the month. (See Appendix C) 

17. City Call Number 

Comment: If the analysis site is located at a city listed in 

the November 1978 Input Data For Solar Systems that site is 
entered into the f-Chart data record. If the analysis site 
Is not a part of the data record, an Interpol atlve routine 
computes the data for any arbitrary site from nearby sites 
where data is available. 

18. Thermal Print Out by Month 
Comment : None 

19. Economic Analysis 

Comment: In general, all runs made for Final Reports specify 

print out of economic analysis. 
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20 . 


Use Optimized Collector Area s 1, Specified Area = 2 


Comment: In general the runs made for Final Reports use 

an optimized collector area. 

21. Solar System Thermal Performance Degradation 
Comment: A value of zero percent Is used. 

22. -46. Economic Parameters 

Comment: The values of the economic parameter were worked 

out between MSFC and IBM for the Final Reports. The source 
of the value is given in the notes on page A-ll. 


Residential 


Item 

Variable Description 

Value 

Units 

Source 

22 

Period of Economic Analysis 

20 

Yrs. 

SAI 1 

23 

Collector Area Dependent System Costs 



MSFC 2 

24 

Constant Solar Costs 



MSFC 2 

25 

Down Payment (% of Original Investment) 

20 

% 

SAI 1 

26 

Annual Interest Rate on Mortgage 

13.5% 

% 

MSFC 2 

27 

Term of Mortgage 

20 

Yrs. 

SAI 1 

28 

Annual Nominal (Market) Discount Rate 

8.5 

% 

SAI 1 

29 

Extra Insur. , Maint. in Year 1 

0.5 

% 

MSFC 2 


(% of Orig. Inv.) 




30 

Annual % Increase In Above Expenses 

10.0 

% 

MSFC 2 

31 

Present Cost of Solar Backup Fuel (BF) 



Actual 

32 

BF Rise: %/Yr. = 1, Sequence of Values = 2 

1 
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Residential (Continued) 


I tern 

Variable Description 

Value 

Units 

Source 

33 

Annual Rate of 6F Rise 





Electricity 

12.5 

% 

MSFC 2 


Oil 

12.5 

% 

MSFC 2 


Natural Gas 

12.5 

% 

MSFC 2 

34 

Present Cost of Conventional Fuel (CF) 



Same as #31^ 

35 

CF Rise: %/Yr. * 1, Sequence of Values - 2 

1 



36 

Annual Rate of CF Rise 





Electricity 

12.5 

% 

MSFC 2 


Oil 

12.5 

% 

MSFC 2 


Natural Gas 

12.5 

% 

MSFC 2 

37 

Economic Print Out by Year = 1, 

2 


Analyst 


Cumulative = 2 



Option 

38 

Effective Federal State Income Tax Rate 





Residential 

30 

% 

SAI 1 


Commercial 

48 

% 

MSFC 2 

39 

True Property Tax Rate Per $ of Original 

0 

% 

SAI 1 


Investment 




40 

Annual % Increase In Property Tax Rate 

NA If 

#39 Is "O' 


41 

Calc. Rt. of Return on Solar Investment? 



Analyst 


Yes = 1 , No = 2 




42 

Resale Value (% of Original INvestment) 

0 


MSFC 2,5 

43 

Income Producing Building, Yes * 1, 



Site 


No = 2 



Dependent 
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Residential (Continued) 


Item 

Variable Description 

Value 

Units 

Source 

44 

Dprc. : Str. In. * 1, Dc. Bal. = 2, 
Sm-yr.-Dgt. » 3, None = 4 

2 


MSFC 2 

45 

If 2, What % of Str. Ln. Dprc. Rt. Is 
Desired 

150 

% 

MSFC 2 

46 

Useful Life for Deprec. Purposes 

20 

Yrs. 

MSFC 2 


NOTES : 

1. Source was Science Applications, Inc. (SAI) Draft Final Report 
on "Comparison of Solar Heat Pump Systems to Conventional 
Methods for Residential Heating, Cooling, and Water Heating," 
April 1979. 

2. These Items were based on judgment and best experience. 

3. The actual current utility rates for the analysis sites 
selected were obtained. (See Appendix D). 

4. The assumption for final report analysis was that the backup 
system actually used for the installation was the same type 
of system that would be used If the solar system was not 
Installed. 

5. The declining balance technique never permits 100% depreciation 
of the asset no matter how long the period. The balance re- 
maining at the end of the system lifetime was treated, for 
accounting purposes, as salvage value. No other salvage 
value was presumed to exist. 
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47. and 48. Economic COPs for Auxiliary Systems 


Comment: These parameters are defined for f -Chart to 
account for economic analysis of solar systems having 
auxiliary backup other than electrical resistance heat. 
The default values of these parameters are as follows: 

Heat Pump Auxiliary COP = 2 

Fossil Fuel Auxiliary COP =0.6 

Electric Resistance COP = 1.0 


The value of the basic COPs are modified, according to the method described 
on page A-2, to account for differences between the fuel used for the 
domestic hot water and the fuel used for space heating. 
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APPENDIX B 

ECONOMIC UNCERTAINTY ANALYSIS 
EQUATIONS 
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APPENDIX B 


ECONOMIC UNCERTAINTY ANALYSIS EQUATIONS 


1. Area dependent Investment costs (C A ) 

alccs ca = -P^ uc A ) 

2. Area Independent Investment costs (C E ) 


ALCCS 


CE 


■P 2 (ac e ) 


3. Ratio of downpayment to Initital investment (D) 

4 lccs d = -(c A A + c E )j pfc g; - 4 j + 


tf(N 


, i, d) [ 


1 - 


fIN.VTJ' 


(AD) 


4. Ratio first year's mlsc. costs to Inlt. Inv. (M) 


ALCCS, 


M 


(C a A + C E ) [ (1 - CT) f (N, g, d) ] (AM) 


5. Ratio first year's assessed value to Inlt. inv. (V) 


ALCCS, 


= - (C a A + C e ) [t ( 1 - D f (N, g, d) ] (AV) 


6. Ratio salvage or resale value to Init. Inv. (6) 


ALCCS, 


■(CflA + C E ) [~ — -v 

A b L o + dr 


(AG) 
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7. Annual market discount rate (d) 


ALCCSd = C FE (L E + L F /C0P F )F(1 - Ct) f(N, e, d) (Ad) 

’(C A A + C E^ | f(N, ] 0? i) Id f ( N ’ °» + 

[d - Ct) M + t (1 - t)vj|g- f(N» g, d) - 

^ " D ) * [f(N,!o, i) Id f ^ N * °* d ^ + 

(* " W 0, Tf)h f(N ’ U d) ] + ( 1 + J jN+1 

* IT" fcT °* d ^ | ^ d ^ 

9 

8. Annual market rate of fuel price increase (e) 

ALCCS e * C fe (L e + L F /C0P F )F (1 - Ct) f(N, e, d) (Ae) 

9. Annual interest rate on mortgage (i) 

iLCCS, = -(C A A + C £ ) |(D-l)(l-t) 

lr f(Hi ') <’ - D > |j - f (w] a; T ) ] 

Iy f (M, 1, d) - t (1 - D) f(N, i, d) 

[' + TWToTT) 2 fr f (M- 0. 1)] J if 
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10. Annual rate of general inflation (g) 


ALCCS g = -(C a A + C e ) ^ (1 - Ct) M + (1 - t) t V 

fg f (N, g, d) (Ag) 

11. Effective income tax rate (t - ) 

ALCCSg = " c FE (L E + L F/ cop F )FCf(N * e * d) ^ 

(C A A + C E ) | (D- 1) t (0M) f(N> u d) 

i_ f WToTT) -t Vf(N, g, d)-C Mf(N, g, d) + 

^f(N, 0, d) | (At) 

12. Property tax rate (t) 

AlCCS t = -(C a A + C E ) (1 - t) Vf(N, g, d) (At) 

13. Cost of electrical energy in the first year (Cp E ) 

alccs cfe = p 1 (l e + l f /cop f )f(ac fe ) 

14. Annual hot water load (L E ) 

ALCCS le = P i C FE F(AL E ) 

15. Annual heating load (l_ F ) 

aLCCS lf = P 1 (C FE /COP F )F(ALp) 
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17. Coefficient of Performance 


ALCCSqqpp = -P] (LpCpp/COPp)F(ACOPp) 

18. Annual load fraction supplied by solar (F) 

ALCCSp = p 1 C FE (L E + L F / cop f^ f) 

NOTE: Three functions used above require definition, as follows: 

f(N, a. b) ■ FTT [ 1 -'(rH) N _ 

i_ f(N , a , b , = [ m. a. b) - rfj (Hi) N ] 

If a ’ = b - a 1 + b ( 1 + b ) “^ N » a ’ 
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APPENDIX C 


MONTHLY AVERAGE WATER 
SUPPLY TEMPERATURES 
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APPENDIX D 


ENERGY COSTS FOR 
ANALYSIS SITES 
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AKRON, OHIO 
GAS 

NOT APPLICABLE 
ELECTRICITY 


0.03854 $/kWh STRAIGHT RATE 

SERVICE CHARGE 6.00$/M0NTH 
ALSO FUEL SURCHARGE 0.0152 $/kWh 
TAX 0 


1000 kWh EFFECTIVE 
RATE* 0.059740 $/kWh 


FUEL OIL 


NOT APPLICABLE 
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ALBUQUERQUE, NM 


GAS (1-505-247-4711) (RESIDENTIAL) 


0-165 THERMS 0.0803/THERM 
165-340 THERMS 0.0826/THERM 
340+ THERMS 0.0966/THERM 

SERVICE CHARGE $1.25 EXAMPLE 

FUEL ADJUSTMENT $0.21 14/THERM 30 THERMS * 0.2114 - 

TAX 4% 


ELECTRICITY (1-505-842-9390) (RESIDENTIAL) 


0-200 kWh 
200-800 kWh 
800+ kWh 
OR 

800 + kWh 


0.05294/ kWh 
0.04794/kWh 
0.03894/kWh 

0.04094/ kWh 


1000 kWh EFFECTIVE 
NOV-MAY RATE = 0.069576 $/kWh 
YEAR-AROUND 

JUN-OCT 


FUEL RATE ADJUSTMENT $0.01 6680/kWh 
SERVICE CHARGE $2.60 
TAX 4.5% 


FUEL OIL 


$0. 999/GAL + 4% TAX 
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FORT WORTH, TEXAS 
GAS 


O-IOOO MCF 
1000-MCF 

SERVICE CHARGE 
TAX 

ELECTRICITY 

0- 25 kWh 
25+ kWh 
FUEL CHARGE 
SALES TAX 


$4.05/MCF 

$2.433/MCF 

0 

0 


$6.00 (MINIMUM) 
$0.0285/ kWh 
$0.008899/ kWh 
4% 


MCF 


1000 kWh EFFECTIVE RATE = $0. 0444/kWh 
FUEL OIL 


NOT USED IN FORT WORTH AREA 



MADISON, WI 


GAS 

0-20 THERMS 
20-50 THERMS 
50+ THERMS 


$0. 28732/THERM 
0.27936/THERM 
0.26892/THERM 


FUEL RATE CHARGE 
ALSO TAX 

SERVICE CHARGE 


$0. 0762/THERM 

0 . 

$2.00/M0NTH 


ELECTRICITY (RESIDENTIAL) 


0- 100 kWh 
100- 500 kWh 
500-1000 kWh 
1000+ kWh 


$0. 0360/kWh 
0.0350/kWh 
0.0320/ kWh 
0.0275/ kWh 


$0. 04167/kWh 


FUEL RATE CHARGE (JAN) $0.00607/ kWh 
ALSO TAX 0 

SERVICE CHARGE $2.00/M0NTH 

1000 kWh EFFECTIVE RATE = 


FUEL OIL 


$0.91 9/GAL 

TAX 0 FOR RESIDENTIAL 4% FOR COMMERCIAL 
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WASHINGTON, DC 
GAS 


$S.OO/MO SERVICE CHARGE 1 THERM = 100,000 Btu 

$0. 3255/THERM + 5 % TAX 

ELECTRICITY (RESIDENTIAL RATES) 

$5.00/M0 SERVICE CHARGE 


NOV - MAY JUNE - OCT 

WINTER RATES SUMMER RATES 


0 - 

600 

kWh 

0.06024 

$/kWh 

0 - 600 

0.06024 

$/ kWh 

600 - 

1500 

kWh 

0.05334 

$/kWh 

600 - 1500 

0.06924 

$/ kWh 

1500 + 


kWh 

0.04289 

$/kWh 

1500 + 

0.26638 

$/kWh 


TAX 16% OF FIRST $15.00 ($2.40 MAX) 

FUEL CHARGE 0.01500 $/kWh (INCLUDED IN ABOVE RATES) 

1000 kWh EFFECTIVE RATE = 0.0675 $/kWh YEAR-ROUND 

FUEL OIL 

$0. 989/GAL + TAX 5% 
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APPENDIX E 


DETERMINATION OF ENERGY 
LOSS (UA) COEFFICIENTS 
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DETERMINATION OF THE UA VALUE OF DETACHED ONE AND TWO FAMILY DWELLINGS 
(Al) AND ALL OTHER RESIDENTIAL BUILDING 3 STORIES OR LESS 

1 . WALLS 


a. Determine the gross area of all exterior walls, including 
windows and doors. (Al 

W 

b. Refer to Figure E-l [8] to obtain combined thermal transmittance 
value (U value) for geographic region. 

OW 

c. Multiply gross wall area by value found in (b) to derive 

u oA for wal1s - 

2. CEILING 

a. Determine total interior surface of ceiling. 

b. For geographic areas where: 

<8000, U QC = 0.05 BTU/H-°F-FT 2 
> 8000, U QC = 0.04 BTU/H-°F-FT 2 

interior ceiling area by value found in (b) to derive 

3. FLOORS 

a. FLOORS OVER UNHEATED SPACES 

(1) Determine the interior floor area (Ap) 

(2) Refer to Figure E-2 to obtain thermal transmittance 
value (Uqp value) in geographic region. 


• HDD 

• HDD 
c. Multiply 
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(3) Multiply interior floor area by value found In (2) to 
derive UppAp for floors. 

b. SLAB ON GRADE FLOORS 

(1) Determine the perimeter of the exposed edge of the 
floor. 

(2) Multiply perimeter length by a factor determined from 

the following table to derive C : L for floor. 

HL r 


t d c hl 

Outdoor Design Heat Loss 

Temperature (°Fl Coefficient (BTU/H-FT) 


-20 to -30 50 

-10 to -20 • . 45 

0 to 10 40 

Above 1 0 35 


(3) Divide the C^Lp product by the difference of the 

outside design temperature (Tp) and the average 
winter building temperature (Tg). 

4. BUILDING UA FACTOR 

The UA factors determined in Steps (1) - (3) are added as follows: 


UA = U 



+ u oc A 


+ U 0F fl F (or C hl L f /(T b - Tp)) 


5. If the UA factor for the building at the actual site is known, computing 

the UA factor as described in Steps (1) - (4) will give a comparison 
value. If this comparison value is less than the given value at the 

actual site, the given value should be used in f-Chart, and the computed 

value for every other analysis site should be increased by the percentage 
difference from the computed value at the actual site. Similarily, if 
the comparison value is greater than the given value for the actual site, 
the given value should be used, and the computed value for every other 
analysis site should be decreased by the percentage difference from the 
computed value at the actual site. 
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U. BTU/H-FT- 


Figure E- 1 


UoWALLS— TYFE“A" BUILDINGS 

TYPE A BUILDINGS SHALL INCLUDE. 

A 1 DETACHED ONE AND TWO FAMILY DWELLINGS 

A 2 ALL OTHER RESIDENTIAL BUILDINGS, THREE 
STORIES OR LESS, INCLUDING BUT NOT LIMITED 
TO: 

MULTI-FAMILY DWELLINGS 
HOTELS AND MOTELS 

ANNUAL CELSIUS HEATING DEGREE DAYS (18 C BASE) 
(IN THOUSANDS) 

12 3 4 5 6 



ANNUAL FAHRENHEIT HEATING DEGREE DAYS (65 F BASE) 

(IN THOUSANDS) 
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BTU/H-FT 1 - 


Flqure F- 2 

Uo VALUES— FLOORS OVER UNHEATED SPACES 

ANNUAL CELSIUS HEATING DEGREE DAYS (18 C BASE) 
(IN THOUSANDS) 

1 2 3 4 5 6 



ANNUAL FAHRENHEIT HEATING DEGREE DAYS (65 F BASE) 

(IN THOUSANDS) 
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APPENDIX F 


ANALYSIS OF AUXILIARY SYSTEM 
OFF-PEAK COOLING PERFORMANCE 
FOR 

SOLARON AKRON 
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ANALYSIS OF AUXILIARY SYSTEM OFF-PEAK COOLING PERFORMANCE 

FOR SOLARON AKRON 
JUNE - SEPTEMBER 1979 


During the 1979 cooling season the total cooling load measured at the 
source was 2.857 million Btu, the sum of the direct and off-peak loads. 
Off-peak cooling furnished 34 percent of the load. The average heat 

pump COP, while operating in the air conditioning cycle, was 2.47 in 
the direct mode and 1.19 in the off-peak mode. 

The energy removed from the off-peak storage tank and transferred to the 
outside environment by the heat pump was 0.554 million Btu during the 
four month season. The energy transferred to the off-peak storage tank 
from the building was 0.983 million Btu. The change in stored energy in 
the off-peak tank was 0.241 million Btu and the first and last off-peak 
storage tank temperature were 49°F and 62°F, which corresponds to the 
energy gain. 

The operating energy required for the transfer of energy from the building 
to the outside environment in the direct mode was 1.007 million Btu. The 
operating energy required in the off-peak mode was 1.922 million Btu, which 
is the sum of the operating energy required to chill water in the off-peak 
storage tank, and to transfer heat from the building to the chilled water. 

The purpose of the off-peak storage system is to take advantage of the 
off-peak electrical rates to pay for storing chilled water which may then 
be used to cool the house with a significatnly lower operating energy paid 

for at peak rates. Performance analysis can be accomplished by comparing 
the amount of energy transferred in the direct mode and in the off-peak 
mode with respect to the operating energy expended in each mode. 

A total of 1.874 million Btu of energy was transferred in the direct mode 
by expending 1.007 million Btu. In the off-peak mode 0.983 million Btu 
of energy was transferred by expending 1.922 million Btu. The efficiency 
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in the off-peak mode was 27.5 percent of the efficiency in the direct 
mode. From the economic point of view 94 percent of the total operating 
energy required in the off-peak mode was paid for at off-peak rates, 
assuming proper operation of the system, and 6 percent was paid for at 

peak rates. To break even, considering the energy transferred, off-peak 
electrical rates would have had to be 27.5 percent of peak electrical rates. 

Overall, the performance of the off-peak system is disappointing. The COP 
of the heat pump operating In the off-peak mode Is 52 percent lower than 

the COP in the direct mode. The pump and blower operate longer to accomplish 
an equivalent energy transfer because of the lower temperature differentials 
achieved in the off-peak mode. The transfer pump more than doubles Its duty 
cycle in transferring energy in the off-peak mode. These energy expenditures 
and system heat gains and losses contribute to make a system that is energy 
inefficient. In addition, there appears to be a low probability of accomplishing 
the purpose of the system, i.e., saving money. 

The data used in this discussion are tabulated in Table F-l. 
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TABLE F-l 


OFF-PEAK COOLING ANALYSIS FOR SOLARON - AKRON 


MONTH 

COOLING 

LOAD- 

DIRECT 

(Mill ion Btu) 

COOLING 

LOAD 

OFF-PEAK 
(Million Btu) 

COP 

DIRECT 

MODE 

COP 

OFF-PEAK 

MODE 

ENERGY 
FROM 
STORAGE 
(Million Btu) 

ENERGY 

TO 

STORAGE 
(Million Btu) 

CHANGE 
IN STORED 
ENERGY 

(Million Btu) 

STORAGE 

TEMP 

FIRST/ LAST 
(°F) 

June 

0.763 

0.575 

2.70 

1.86 

0.457 

0.575 

0.215 

49/69 

July 

0.159 

0 


- 

-0.546 

0 

-0.059 

84/93 

August 

0.646 

0.240 


0.79 

0.382 

0.240 

0.003 

53/57 

Sept. 

0.306 

0.168 

1 

0.92 

0.261 

0.168 

0.082 

54/62 

imonii 

1.874 

0.983 


- 

0.554 

0.983 

0.241 

- 

Average 

0.469 

0.246 

2.47 

1.19 

0.139 

0.246 

0.060 

- 


MONTH 

OPERATING 
ENERGY 
DIRECT MODE 
(Million Btu) 

DIRECT OPERATING 
ENERGY-OFF 
PEAK MODE 
(Million Btu) 

OFF-PEAK 

OPERATING 

ENERGY 

(Million Btu) 

TOTAL OFF-PEAK 

OPERATING 

ENERGY 

(Million Btu) 

— 

mm 

0.065 


0.383 

I 

■SQM 

0 

1 

0.322 

■ 

0.352 

0.028 


0.761 


0.159 

0.022 

0.434 

0.456 

Total 

1.007 

0.115 

1.807 

1.922 

Average 

0.252 

0.029 

0.452 

0.481 
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